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Abstract 

The asymptotic properties of self-similar spherically symmetric perfect 
fluid solutions with equation of state p = afj, (— 1 < a < 1) are described. 
We prove that for large and small values of the similarity variable, z = r/t, 
all such solutions must have an asymptotic power-law form. They are asso- 
ciated either with an exact power-law solution, in which case the a > ones 
are asymptotically Friedmann, asymptotically Kantowski-Sachs or asymp- 
totically static, or with an approximate power-law solution, in which case 
they are asymptotically quasi-static for a > or asymptotically Minkowski 
for a > 1/5. We also show that there are solutions whose asymptotic be- 
haviour is associated with finite values of z and which depend upon powers 
of Inz. These correspond either to a second family of asymptotic Minkowski 
solutions for a > 1/5 or to solutions that are asymptotic to a central sin- 
gularity for a > . The asymptotic form of the solutions is given in all 
cases, together with the number of associated parameters. This forms the 
basis for a complete classification of all a > self-similar solutions. There 
are some other asymptotic power-law solutions associated with negative a, 
but the physical significance of these is unclear. 



1 Introduction 

Spherically symmetric solutions to Einstein's equations admitting a self-similarity 
of the first kind are characterized by the fact that the spacetime possesses a ho- 
mothetic Killing vector. This means that they can be put into a form in which 
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every dimensionless variable is a function of some dimensionless combination of 
the cosmic time coordinate t and the comoving radial coordinate r. Such solu- 
tions have attracted considerable attention in recent decades. This is partly for 
mathematical reasons, since the governing equations reduce to ordinary differen- 
tial equations from partial differential equations. In particular, in the case of a 
perfect fluid, the only barotropic equation of state compatible with the similarity 
assumption is one of the form p = a/x, where a is a constant, and the equations 
then simplify still further. However, it is also for physical reasons, since such 
solutions play a crucial role in many cosmological and astrophysical contexts and 
are important in the context of the "critical" phenomena discovered in recent 
gravitational collapse calculations (Evans and Coleman 1994, Carr ct al. 1999, 
Neilsen and Choptuik 1998). Indeed the "similarity hypothesis" proposes that 
spherically symmetric solutions may naturally evolve to similarity form in a wide 
variety of situations. The arguments for this have been reviewed by Carr & Coley 
(1999a). 

The possibility that self-similar models may be singled out in this way from 
more general spherically symmetric solutions means that it is essential to study 
the full family of such solutions. Two main approaches have been used in such 
studies. The first uses the "comoving" approach, pioneered by Cahill & Taub 
(1971), in which the coordinates are adapted to the fluid 4-velocity vector and 
the solutions are parametrized by the similarity variable z = r/t. At a given 
value of t, this specifies the spatial profile of the various quantities. At a given 
value of r (i.e. for a given fluid element), it specifles their time evolution. The 
second uses the "homothetic" approach, introduced by Bogoyavlensky (1977). In 
this, the coordinates are adapted to the homothetic vector and solutions can be 
treated as orbits in a compactified 3-dimcnsional phase space. 

The purpose of this paper is to delineate all possible asymptotic behaviours of 
spherically symmetric p = similarity solutions using the comoving approach. 
This provides the first step in our "complete" classification of all such solutions 
with « > in another paper (Carr & Coley 1999b; CC). The classification is 
complete, subject to the requirement that the solutions are physical everywhere 
and do not contain shocks. However, our claim for completeness in that paper 
is based on the assumption that all solutions must have an asymptotic form 
in which any dimensionless quantity has a power-law dependence on z at large 
and small values of z or a power-law dependence on Inz at finite z (whenever 
this corresponds to zero or infinite distance from the origin). We do not prove 
this assumption in CC but take it as our starting point. We then find that the 
behaviour of solutions in the asymptotic limits can only take one of a few simple 
forms for a > 0: asymptotically Priedmann, asymptotically Kantowski-Sachs, 
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asymptotically Minkowski (for a > 1/5), asymptotically singular or what we 
term asymptotically "quasi-static" . 

In this paper we also obtain the asymptotic behaviours but in a different way 
and, in the process, we provide a rigorous proof of the "power-law" assumption 
in the physically important positive a case. This paper therefore crucially com- 
plements the classification paper. Indeed wc sometimes refer to equations in CC 
below. Both papers end up with the same asymptotic behaviours for a > but 
the analysis in this paper is more rigorous and more comprehensive. In addition, 
this paper also covers the negative a case. However, wc do not consider the phys- 
ical significance of any of the solutions in this paper. The positive a solutions are 
discussed in more detail elsewhere (CC, Carr et al. 1999). The physical signifi- 
cance of the negative a solutions are less clear; some applications are discussed in 
Carr & Coley (1999a) and we note that for —1 < a < —1/3 the asymptotically 
Friedmann (AF) solutions (for example - see Table 1) are inflationary (Ohve, 
1990) 

It should be stressed that Goliath ct al. (1998a & b; GNU) have also implicitly 
provided a classification of spherically symmetric self-similar solutions with posi- 
tive a using the homothetic approach. However, the nature of their classification 
is rather different in that they emphasize the equilibrium points of the associ- 
ated dynamical system. These always lie on the boundary of the 3-dimensional 
phase space, whereas the condition z ^ oo emphasized by CC usually corre- 
sponds to the interior of the phase space and so does not play a crucial role in 
the GNU analysis (for example, there are no equilibrium points corresponding to 
CC's asymptotically quasi-static solutions which occur as \z\ — > oo). However, 
some of these equilibrium points do have a straightforward connection to CC's 
analysis. Thus what GNU term the "F point" corresponds to solutions which are 
asymptotically Friedmann as 2; — > oo, their "C° point" corresponds to solutions 
which are asymptotically Friedmann as 2; — > (or regular at the origin), their 
"T point" corresponds to the static solution, and their "K point" corresponds to 
(asymptotically Kasner) solutions with a central singularity (occuring at a finite 
value of z). The two approaches are therefore complementary and the value of 
using a combined approach has been emphasized by Carr et al. (1999). 

The outline of this paper is the following. In Section 2 we introduce the 
general spherically symmetric similarity equations for a perfect fluid and highlight 
the features necessary for our analysis. In Section 3 we provide the analysis 
itself; this comprises an exhaustive case-by-case consideration of the different 
possible behaviours of various key functions. In Section 4 we draw some general 
conclusions and point out some more subtle issues which are not covered exphctly 
in Section 3. 
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2 Spherically Symmetric Similarity Solutions 



In the spherically symmetric situation one can introduce a time coordinate t 

such that surfaces of constant t are orthogonal to fluid flow lines and comoving 
coordinates (r, 6, 0) which are constant along each flow line. The metric can then 
be written in the form 

ds^ = e^" df - e^^ dr^ - dfl^, d^f = dO^ + sin" 6 d(f)^ (2.1) 

where u, A and R are functions of r and t. For a perfect fluid the Einstein 
equations are 

G^"" ^8n[{fx + p)U''U'' -pgi"'] (2.2) 

where //(r, i) is the energy density, p{r,t) the pressure, — (e"*^, 0,0,0) is the 
comoving fluid 4- velocity, and we choose units in which c — G — 1. The equations 
have a first integral 



m(r, t) — ^R 



(2.3) 



and this can be interpreted as the mass within comoving radius r at time t: 

dR 

m(r,t)=47r jiR^—dr'. (2.4) 
Jo or 

Unless p = 0, this quantity decreases with increasing t because of the work done 
by the pressure. 

A self-similar solution is one in which the spacetime admits a homothetic 
Killing vector ^ such that 

^^^■,u + ^u;ii — '^9^lu■ (2-5) 

This means that the solution is unchanged by a transformation of the form t — > at, 
r ^ ar for any constant a. Solutions of this type were first investigated by Cahill 
& Taub (1971), who showed that by a suitable coordinate transformation they 
can be put into a form in which all dimensionless quantities such as A, 

S=-^ M=5, P = pR^, W = iJiR^ (2.6) 
r R 

are functions only of the dimensionless variable z = r/t. Then we have 

^ = ^ = (2 7) 

dt r dz^ dr r dz^ 
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so the field equations reduce to a set of ordinary differential equations in z. 
Another important quantity is the function 



V{z) = e^-'^z, (2.8) 

which represents the velocity of the surfaces of constant z relative to the fluid. 
These surfaces have the equation r — zt and therefore represent a family of 
spheres moving through the fluid. The spheres contract relative to the fluid for 
z < and expand for z > 0. Special signiflcance is attached to values of z for 
which |y| = 1 and M = 1/2. The first corresponds to a Cauchy horizon (either 
a black hole event horizon or a cosmological particle horizon), the second to a 
black hole or cosmological apparent horizon. 

The only barotropic equation of state compatible with the similarity ansatz 
is one of the form p = afi (— 1 < a < 1). It is convenient to introduce a 
dimensionless function x{z) defined by 

x{z) = (47r/ir') (2.9) 

The conservation equations T^'-'.^^ = can then be integrated to give 

= (2.10) 

A -l/ag2 (2.11) 



e = 7x 



where /? and 7 are integration constants. The remaining field equations reduce 
to a set of ordinary differential equations in x and S: 



l + ay (5'^ S X \ a 

M = ^2^-(i+")/« 



1 + + 



(2.13) 
(2.14) 



2 

M = - + ^a;-2;22(l-a)/(l+a)^2 _ 1^2^-(2/a)^6 + ^ , (2.15) 

2 2/^ 2 \ S J 

where the velocity function is given by 

V= (/37)-ia;(^-")/"5-'^(i-")/(^+") (2.16) 
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and an overdot denotes zd/dz. 

We can envisage how these equations generate solutions by working in the 
3-dimensional (x, S, S) space (Carr and Yahil 1990). At any point in this space, 
for a fixed value of a, eqns (2.14) and (2.15) give the value of z; eqn (2.13) then 
gives the value of x unless \V\ = ^Jol and eqn (2.12) gives the value of S. Thus 
the equations generate a vector field (i, S*, S) and this specifics an integral curve 
at each point of the 3-dimcnsional space. Each curve is parametrized by z and 
represents one particular similarity solution. This shows that, for a given equation 
of state parameter a. there is a 2-paramctcr family of spherically symmetric 
similarity solutions. Special significance is attached to points where \y\ = ^/a, 
which specifies a 2-dimensional surface in {x, S, S) space. This corresponds to 
a sonic point, so there can be a discontinuity in the pressure gradient (i.e. the 
value of x) here. The requirement that solutions be regular at the sonic point 
(in the sense that they can be extended beyond there) places severe restrictions 
on the nature of the solutions (Bogoyavlensky 1977, Bicknell & Henriksen 1978, 
Carr & Yahil 1990, Ori & Piran 1990) but these restrictions are not relevant to 
the present asymptotic analysis. 



3 Asymptotic Classification 

In this section we will provide a complete analysis of spherically symmetric sim- 
ilarity solutions in the limits i? — > and R ^ oo. (We will assume z > 
throughout since the z < solutions are just the time reverse of the z > ones.) 
This may correspond to z ^ 0, z ^ oo or z ^ z^ (finite) since the correspondence 
between the comoving radial coordinate r and the physical distance R may not 
be straightforward. In particular, r — > need not imply i? — > (eg. some dust 
solutions), R ^ need not imply r — (eg. black hole singularities), r — > oo 
need not imply i? — > oo (eg. the Kantowski-Sachs solution) and R ^ oo need not 
imply r — > oo (eg. asymptotically Minkowski solutions). We categorize solutions 
according to whether y — > 0, F — > oo or y ^ K (finite) in each of the z limits. 
A key role is played by eqn (2.13), which can be written as 

X _ 2V'^a S lo? 1 a(l + g) ^-S^^"")/" 

One has different possible behaviours according to whether ^ — > 0, ^ — > oo or 

^ — (finite). In principle we need to analyse 27 cases, corresponding to the 
three limiting behaviours of each of 2;, V and ^. However, in practice, many of 
these do not lead to a consistent solution. In the discussion below, we usually 
only derive the self-consistent ones. 
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Our procedure is as follows: for each limiting value of z, V and ^, we use eqn 
(3.1) to eliminate x/x in eqn (2.12) and thereby obtain a second order ODE for S. 
Not every solution of this ODE may be consistent with the other equations, so we 
first confirm that ^ has the correct limiting value and then check that the integral 
condition associated with eqns (2.14, 2.15) is satisfied. This procedure not only 
identifies the possible asymptotic behaviours but also indicates the number of 
free parameters in each case. The results are summarized in Table 1. 





z^O 


z ^ z^ 


z —>■ oo 




AKS (0/1) (1 > a > -1/3) 

AF (1) (1 > a > -1/3) 
ES (0) (1 > a > 0) 
AX (1) (-0.17 <a< -1/7) 


No soln 


AF (1) (-1 < a < -1/3) 
AKS (1) (-K a < -1/3 
AX (1) (1 > a > -1/7) 




No soln 


AM (2) (1 > a > 1/5) 
AKS (1) (1 > a > 0) 


AM (1) (1 > a > 1/5) 
AM (1) (-K a < 0) 


V^oo 


AKS (1) (-K a < -1/3) 
AF (1) {-l<a< -1/3) 


AK (2) (1 > a > -1/3) 


AQS (2) (1 > a > 0) 
AKS (0/1) (1 > a > -1/3 
AF (1) (1 > a > -1/3) 
AY (1) (-K a < 0) 



Table 1. This summarizes the asymptotic (A) behaviour of spherically symmetric 
perfect fluid similarity solutions for different ranges of a. The number of arbitrary 
parameters for each solution is given in parentheses. KS, F, M, K, ES and 
QS denote Kantow ski- Sachs, Friedmann, Minkowski, Kasner, Exact Static and 
Quasi-Static respectively. X and Y correspond to the two new solutions in which 
either V or M is negative. The AKS solution is described by one parameter 
for q; > and —1 < a < —1/3 but there is only the exact KS solution for 
— 1/3 < a < V and M are positive only for — 1 < a < —1/3. 



In the following analysis a differential equation of the form 

S + aS + b^^Q, (3.2) 

where a and b are constants, will frequently arise. When a 7^ 0, we can divide 
through by S and integrate this equation to obtain 

SS^ ~ z-". (3.3) 
If 6 7^ —1, we can then integrate a second time to obtain the exact solution 
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5 = (ci + C2^-")Va+^), (3.4) 

where Ci and C2 are integration constants. The dominant asymptotic solution 
is then determined by the sign of a. Note that the first order perturbation to 
the dominant solution (which is not derived in the present analysis) will often be 
larger than the subdominant solution. 

In this case, eqn (3.1) becomes 

- = -T^-«. i^-'-^S-V-'". (3.5) 
If ^ — > 0, this imphes x ~ ^-'^a/i^+a) ^ gQ (2.12) gives 

sU'-^)s+'4 + ^^0. (3.6) 

\l + a J S 1 + a ^ ^ 

Defining a{z) = zS{z), this becomes 

1±1^^)*+H^!=0. (3.7) 

This differential equation is of the form of (3.2) and, from eqn (3.4), we therefore 
obtain the exact solution 

For ci = 0, we obtain the exact power-law solution 

where the constants A and B are constrained by eqns (2.14, 2.15) and can be 
taken to be 1 with a suitable choice of the scaling constants /? and 7 in eqns 
(2.10, 2.11) [see CC eqn (3.10)]. This corresponds to the k=0 Priedmann model. 
Note that F — > and ^ — > (as assumed) providing a > —1/3. As indicated 
in Tabic (1) and discussed in CC, there is also a 1-parameter family of solutions 
asymptotic to this. This involves modifying the constants A and B, as well as 
adding higher order terms in z. However, demonstrating this depends on a higher 
order analysis than presented here. 

For C2 = 0, one obtains another power-law solution 

S = Az-\ a; = (3.10) 



8 



which corresponds to a Kantowski- Sachs (KS) model. Eqns (2.14, 2.15) determine 
the constants A and B in terms of a once (3 and 7 are specified [see CC eqn (3.18)]. 
However, as explained in CC, one requires /3 and 7 to be imaginary for a > 0, 
which means that V is negative (i.e. the solution is tachyonic). Furthermore, 
eqn (2.14) implies that M is negative for a > —1/3, so the solution given by eqn 
(3.10) may be unphysical. Nevertheless, V ^ and ^ ^ providing a > —1/3, 
so this solution is still formally valid. As indicated in Table (1) and discussed 
in CC, one can perturb this solution to obtain a 1-parameter family of solutions 
which asymptote towards it. This involves modifying the constants A and B, as 
well as adding higher order terms. If both constants Ci and C2 in eqn (3.8) are 
non-zero, the Ci term in eqn (3.8) dominates as z but the other term should 
not be interpreted as giving the perturbation to KS. 

If ^ — > ^* (constant), then x ~ and eqn (2.12) becomes 



y 1 — Q!^ J \ 1 — a J S 

This is again of the form (3.2) and hence has the exact solution 

S=[c^ + c2;,(i+6"+"^)/(i-"')]("-i)/(i+7«). (3.12) 

If C2 = 0, this leads to the exact static solution with (cf. Misner & Zapolsky 
1964) 

X = Xo{a), S = 5*0(0;). (3.13) 

This solution is specified uniquely for a given equation of state since eqn (3.5) 
and eqns (2.14, 2.15) give two independent relationships between Xo and So- [See 
CC eqn (3.30) for the explicit form of eqn (3.13).] However, these values are real 
and the mass is positive only for ct > 0. 
If Ci = 0, then 

g ^ ^^-(l+6a+a2)/(l+a)(l+7a)^ ^ ^ ^-4a{l+6a+a^)/{l-a'^)il+7a) ^ (3-14) 

where the constants A and B are related by eqns (2.14, 2.15) and by eqn (3.5). 
From eqns (2.14) and (2.16) one also has 

V ~ ^-(3«+l)V(l+")(l+7a)^ M ~ ^2(3a+l)(l+6a+a2)/(l-a2)(l+7a)^ ^3_;^5^ 

SO this is only consistent with the assumption that V^Oas^;— *Oif— 1<q;< 
-1/7. However, eqns (3.5) and (3.14) require -2q;(1 + 3af/{l - a'^){l + 7a) = 
(1 + 0)7-2 A-'^E-^^-")/", and this is only satisfied for -I < a < -1/7 if 7 (and 
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hence f3) are imaginary. In this case, M is positive but V is negative (as in the 
KS case), which might be regarded as unphysicaL However, one can show that 
eqns (2.14, 2,15) are still satisfied, and indeed give the same relationship between 
A and B as above, providing — 1/7 > a > 2^/2 — 3 = —0.17. Since this does not 
seem to correspond to any well-known solution, we describe it as asymptotically 
"X" in Table (1). If both constants Ci and C2 in eqn (3.12) are non-zero, then 
the Ci term dominates as z ^ for a > 2^2 - 3. Although one might regard 
this as an asymptotically static solution, it is not equivalent to a perturbed static 
solution. 

If ^ ^ oo, it can be shown that there is no consistent solution. 



In this case, eqn (3.1) implies 

- = 2«l + I^T^I - e = «(1 -h a)/52x("-^)/"^-^(^-)/(^+"), (3.16) 

where we have used eqn (2.16) to eliminate V in the expression for ^ The term in 
curly brackets is included to allow for the possibility that x/x — > and S/S ^ 0. 
If ^ ^0, this imphes x ~ S"^" and so eqn (2.12) becomes 

Eqn (3.2) implies that the solution of this equation is given by 

5= [ci + C2^("-l)/(l+")]V(l-2-). (3.18) 

The case C2 = gives the the exact static model, in which case Xo and So have 
the values indicated by eqn (3.13) and one needs the term in curly brackets in 
eqn (3.16). The case ci = gives another 1-paramctcr solution, 

g ^ j^^{a-l)/{l+a){l-2a) ^ ^ ^ ^^2a(a-l)/(l+a)(l-2a) ^ (g -^g) 

where A and B are related by eqns (2.14, 2.15). However, this leads to 

\/ rs. ^(l-")/(l+")(l-2«)^ ^ ^ ^2a(l-a)/(l+a)(l-2a)^ (3.20) 

SO y ^ oo and ^ — > as z — >■ cxd only for a < 0. In this case, eqns (2.14, 2.15) 
can be satisfied only if f3 and 7 arc imaginary, so that V is negative, but M is 
still positive. Since this does not seem to correspond to any well-known solution, 
we describe it as asymptotically "Y" in Table (1). 
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For ci 7^ and C2 7^ 0, eqn (3.18) leads to a solution 

S^A + , x^C + . (3.21) 

In this case, the term in curly brackets in eqn (3.16) can still be neglected since 
x/x and S/S go to zero more slowly than V~^. There are two relationships 
between the four integration constants A, B, C, D, eqn (3.21) requiring D — 2aB 
and cqns (2.14, 2.15) giving another relationship, so these solutions are described 
by two independent parameters. Note that S tends to a constant as 2; 00 but 
this constant is different from the value 5*0 given by eqn (3.13), so these solutions 
are described as asymptoticahy "quasi-static". Indeed eqns (3.21) corresponds to 
the full family of such solutions discussed by CC and also (implicitly) by Ori & 
Piran (1990) and Foglizzo & Henriksen (1992). 

If ,^ ^ (constant), then x ~ ^-'ia/{i+a) (3.18) again applies. For 

a > —1/3, this leads to cither the Fricdmann solution, given by eqn (3.9), or the 
KS solution, given by eqn (3.10). Although the general solution with ci 7^ and 
C2 7^ can be written as 

S = ^^-2/3(l+«) [1 + 5^-(l+3a)/(l+a)] ^^ ^2) 

and might therefore be regarded as being asymptotic to the Priedmann solution, 
the second term is not equivalent to the first order perturbation of the Fried- 
mann solution. CC show that this pertubation goes like 2;~^(^+3")/3(^+") and this 
dominates the second term in eqn (3.22) as z —>■ oo. 
If ^ — > oo, there is no consistent solution. 



• y— ^ooasz— >0 

In this case, eqn (3.16) apphes. If ^ — > 0, one again obtains eqn (3.18). 
However, the static solution with C2 = has K — > for all a, while eqn (3.20) 
implies that no solution with Ci = has both V ^ oo and ,^ ^ 0. If 
gets both a Friedmann solution of the form (3.9) and a KS solution of the form 
(3.10) with y — > oo and ^ — > providing —l<a< —1/3. Furthermore, unhke 
the a > case, the KS solution is now physical in that M and V are positive. If 
^ — > oo, there is no consistent solution. 

• K— >Oasz— >oo 

In this case, eqn (3.5) applies. If ^ 0, one again obtains eqn (3.8) and one 
gets both a Friedmann solution and a (physical) KS solution with V ^ and 
,^ ^ providing — l<ci;<— 1/3. If^^^*, one gets a solution of the form 
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(3.12). However, the static solution with C2 = has V — >■ 00 for all a, while eqn 
(3.5) cannot be satisfied for ci = because 

X 2a _ 2a{l + 3af 

X l + a (l-a2)(i + 7«) ■ 
and therefore has the same sign as ^. If ^ ^ 00, there is no consistent solution. 

• V ^ 00 a,s z ^ 

In this case, eqn (3.1) implies 

^ = 2af - i = a{l + a)/32a;("-i)/";^;2(i-a)/(i+a) ^ ^3 24) 

X o 

(cf. eqn (3.16) without the term in curly brackets). If ^ — > 00, then x/x, S/S 
and ^ must all diverge at the same rate (i.e. there is no self-consistent solution 
with x/x « S/S or x/x >> S/S). Hence, 

- ^ (3.25) 

X 

which integrates to 

x = (di + dsM"^^'""^- (3-26) 
To ensure ^ — > cxd as 2; ^ 2;^, one requires a; for o; > or x — > 00 for a < 0. 
In both cases, one needs di = —d2lnz^ and hence 

x^[ln{z/z.)]''/^'-''\ (3.27) 

Defining 

L = ln{z/z^), (3.28) 
it is easy to show that the only consistent solution is then 

5 = ^L2/3(i-°), a; = ^^"/(i-'^) (3.29) 

where A and B are constants related by eqns (2.14, 2.15) and eqn (3.24). Thus 
the scale factor goes to zero, while the density diverges for a > and goes to 
zero for a < 0. However, we also have 

V ~ L-(l+3a)/3(l-a) (3 30) 

and this diverges as 2; — > (as required) only for a > —1/3. Eqn (2.14) implies 

M~L-2/3(i+"), (3.31) 
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so M — * oo at z = 2;* and MS constant. Eqn (2.14) also implies that M 
is positive providing the sign of S/S is positive, so this requires that z —>■ z^, 
from above. Note that eqn (3.24) yields the same relation between A, B and z^ as 
eqns (2.14, 2.15), so these solutions are described by two independent parameters. 
Eqns (2.10) and (2.11) imply that the metric components are 

~ ^0, ~ l-a/3{l-a) ^ ^ ^2/3(l-a) ^ q_ (3 32) 

corresponding to a singularity of infinite density (cf. Schwarzschild). 

If ^ — > or (constant), it can be shown that there is no consistent solution. 



• y — > as 2; — > 

In this case, eqn (3.1) implies 

X ^ 2V^a S a(l + «)^ -2(l-a)/(l+a)^(a-l)/a 

X V-aS l + aV,^-a ^' ^ " V,^ - a 

(3.33) 

where V rather than V^, appears in the first term on the right because the product 
{V^ — V^){S/S) may tend to a finite hmit. If ^ — > 0, this implies 

X ~ 52y>/(K^-a) (3 34) 

providing x/x and 5"/ 5" diverge. However, unless V /V — > oo, eqn (2.16) also 
imphes x ~ so we require — 1. Eqn (2.12) then becomes 

/l _4q;_Q,2\ . g2 

S + " \s--^^^0. (3.35) 

y 1 — J 1 — a b 

This is of the form of eqn (3.2), so the exact solution for o; 7^ 1/5 and a 7^ \/5 — 2 
is given by 

S^[ci + rf(^/^^){"^+4a-l)/{l-«2)]{l-«)/{l-5a)_ (333) 

In order to have ^ — > as z — > 2;*, eqn (3.33) implies that one needs S to go to 
infinity, whatever the sign of a. One therefore requires d — —Ci and a > 1/5. 
This excludes a < and also the special cases o; = 1/5 and a — y/b — 2. Noting 
that to first order 

we obtain the asymptotic form 

5 = AL(^-°)/(^-^"), a; = 5L2"/(i-5°) (3.38) 
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where A and B are constants. Thus the scale factor diverges and the density goes 
to zero. The condition = 1 gives a relationship between the constants A, B 
and z^, so these solutions are described by two parameters. 
Eqn (2.14) implies that 



M ~ /,(l-")/{5a-l) 



(3.39) 



so these solutions have zero mass ai z = z^, and MS tends to a constant. Fur- 
thermore the sign of S/S implies that z must approach z^ from helow in order 
that the mass be positive. Eqn (2.15) can be written as 



M 




(3.40) 



Since x'^'^'S^ - ^2(3a-i)/(5a-i) gQgg infinity for a < 1/3 and zero for a > 1/3, 
one requires the term in curly brackets to go to to zero and infinity, respectively, 
in these two cases. However, the last term in eqn (3.40) also scales as 



^(a-l)/(5a-l) 



2-£ 

S 



(3.41) 



so we need the term in square brackets to go to zero as L^^ a)/{5a i) both cases. 
We therefore need 

^iV - 1) - 2. (3.42) 



Using eqns (2.16) and (3.33), together with the relation 

V 1 



one then obtains 



— a 1 — a 

V _ 

V ~ 



1 



a 



1 — a 
1 — 5a 



iy^ - 1) 



1 



a 



(3.43) 



(3.44) 



Condition (3.41) also determines the second order terms in the expressions for x 
and S but does not impose any further relationship between A, B and z^. Note 
that the metric can be written as 



(3.45) 



This resembles the open Friedmann model and is related to Minkowski space- 
time by a time-transformation (as in the Milne model). These solutions can 
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therefore be regarded as asymptotically flat, or, more precisely, as aymptotically 
Schwarzschild since the mass (~ MS) tends to a non-zero constant. 

If ^ — > oo or constant, it can be shown that there is no consistent solution. 



In this case, we require y — > and eqn (2.16) then implies 

S 1/1 — Q;\i; 1/1— q; 



(3.46) 



S 2 \ a J X 2 \1 + a. 
The exact solution of this differential equation is 

^2 ^ SozT^x'-^ (3.47) 

which can then be combined with eqn (3.1) to give 

x-(K^-l)(l + a) ^' V^-l ^'-^'^ 

where ^' differs slightly from ^. If ^' — > 0, this implies 

x^Bz- ^ = l^il^^l±^ (349) 
and eqn (2.16) then gives 

where A and B are integration constants. For a > 0, these solutions are consistent 
with the condition ^' — > as 2; — > cxd providing x ^ 00 and, from eqn (3.48), this 
requires > 1. This in turn imphes S ^ 00. Eqn (2.14) gives 

^ ^ {V^ -a){l + a) ^_^y2(i_^)+i+3^l/(y^2_i^(i^^^ 

2a{V^ - 1) ^ ■ ^ 

and so the mass is positive and tends to zero for > 1. On the other hand, 
eqns (2.15), (3.48) and (3.49) imply 

M - i ~ 1 1^1^ (V^ (3.52) 
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If the exponent of z in this expression is positive, M can go to zero only if the 
term in curly brackets does and this requires SjS ~ Equating this 

with the value oi SjS imphed by eqn (3.48) gives a quadratic equation for V^: 

(1 - a)Kf - 2a(l + a)K - a(l + 3a) = (3.53) 

with the solution 



a{\ + a) + \ aid? — + 3q; + 1) 
K = ^ , -. (3.54) 

As a. decreases from 1 to 0, K decreases from oo to and it exceeds 1 (as required) 
only for a > 1/5. Note that eqn (3.52) implies that the exponent of z in eqn 
(3.51) is indeed positive (as assumed). 

The value of is also well-defined for — 1 > o; > —1/3. However, if one seeks 
solutions of this kind, then the condition ^' — > as ^ — > oo implies a; — > and eqn 
(3.48) would then give 1 > > —2a /{I — a). The lower limit is incompatible 
with eqn (3.53), so there are no a < solutions. One can also formally obtain 
a negative value of V* by taking the negative square root in eqn (3.53). This is 
relevant if one is considering asymptotically KS solutions with a > 0. 

Eqns (3.48), (3.49) and (3.53) impose a relationship between the constants A 
and B and so these solutions are described by one-parameter and they all tend 
to the same value of V for a given a. The requirement that the term in curly 
brackets in eqn (3.52) goes to zero also determines the second order terms in the 
expansions for x and S. These solutions, like the previous ones, are asymptotically 
flat: the metric components can be expressed as 

e-^^Km^-i), e^^zV(^'-^), (355) 

and the metric can be reduced to the Minkowski form with a suitable change of 
coordinates. 

If ^ — > oo or constant, there is no consistent solution. 



• V ^ V^. as z ^ 

In this case, eqns (3.46) - (3.53) still apply but, for a > 0, the condition 
^' ^ requires < 1. Eqn (3.50) then implies that the mass is negative unless 

< 1/y/a. However, eqn (3.53) requires that one always has V* > 1/y/a, thus 
excluding this case. For a < 0, the condition ^' — > requires either > 1 or 

< —2a/{l — a). The first possibility is excluded because the mass is negative 
but, the second possibility is allowed providing —1 < a < —1/3, so that V* is 
defined. 



16 



4 Discussion 



In this paper we have focussed on estabhshing the asymptotic forms of spherically 
symmetric perfect fluid similarity solutions. The results are summarized in Table 
1. Only those cases that give rise to self-consistent and physically reasonable 
solutions have been described explicitly (none of the remaining cases give rise to 
self-consistent solutions). Although some of the details of the calculations have 
been omitted, the analysis presented can be regarded as a rigorous justification for 
the underlying assumption of our complete classification of positive a solutions 
in Carr & Coley (1999b). Although our description of these solutions has been 
rather brief, their physical properties have been further studied by Goliath et al. 
(1998a, 1998b) and Carr et al. (1999). Our analysis has covered all equations of 
state with |q;| < 1 but the stiff fluid case (a = 1) has not been exphcitly included 
(although many of the results carry through in this case); this is because different 
mathematical techniques are more appropriate in this case and the analysis will 
be presented elsewhere. 

It should be stressed that our analysis is not completely rigorous since we have 
implicitly assumed that certain regularity conditions are satisfied. For example, in 
Section 3 we always presupposed that V has a well-defined asymptotic limit. Also 
our analysis has assumed certain properties for the derivatives of small quantities. 
However, an asymptotic relationship of the form X ^ Y does not necessarily 
imply X ^ Y asymptotically. In order to "prove" the results of Section 3 we 
must show that this type of behaviour is not possible and this must be done on 
a case-by-case basis. For example, suppose an asymptotic solution of the form 
X ~ 2;Ms deduced. In principle 

xr^z\\+e) (4.1) 

where 

- = / + E=-^. (4.2) 

However, the condition e ^ 1 does not necessarily imply E <^\. To prove the 
result, we must keep e and e in all of the equations and show that the only self- 
consistent solutions must have £ -C 1. Suffice it to say that all cases have been 
checked and the only self-consistent solutions are those given in this paper. In 
fact, studying solutions of the form given by eqn (4.2) is precisely what we do in 
CC in deriving solutions which are asymptotic to the Friedmann and Kantowski- 
Sachs models (although in practice this calculation was divided up into different 
separate cases than in that of section 3-and this served as a further check of the 
validity of our results) . 
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The results of Section 3, and the existence of well-defined limits and sufficient 
regularity, also follow from the mathematically more rigorous dynamical systems 
analyses of Bogoyavlensky (1985) and Gohath et al. (1998a, 1998b), at least in 
the positive a case. In these analyses, monotone functions and Dulac functions 
were found to exist, thereby prohibiting periodic orbits and limit cycles in the 
corresponding phase spaces. This rules out the possibility of kinds of asymptotic 
behaviour different from those discussed in this paper. Also, many of the results 
apply to the case o; = and are consistent with the results of the analytic study of 
the dust models (Carr 1999); indeed, this serves as a further check of the validity 
of our results here. 

Finally, we should emphasis that our analysis has demonstrated that all self- 
similar spherically symmetric perfect fluid solutions are asymptotically of power- 
law form for z — > oo or 2; — > and of log-power-law form for 2; — > 2;*. This 
provides the basis of the work presented in our classification paper. 
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